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Abstract 


The  diffraction  of  a  plane  wave  "by  an  infinite  grating 
of  identical  circ-'olar  cylinders  is  considered.  It  is  shown  that 
the  various  scattered  cylircxical  waves  combine  to  form  c  finite 
number  of  plane  waves  and  an  infinite  miniber  of  "  STirface  waves" 
propagated  along  the  grating. 
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I.  Introduction 

We  consider  the  diffraction  of  a  plane  wave  ty  an  infinite  grating 
of  identical  circular  cylinders.  The  totsil  scattered  field  can  "be  represented 
as  a  series  which  is  the  sum  of  waves  scattered  "by  the  individual  cylinders, 
We  shall  not  he  concerned  with  the  coefficients  in  the  representations  of 
these  waves,  hut  rather  consider  only  the  summation  of  the  waves  over  the 
individual  cylinders  of  the  grating.  We  show  that  the  scattered  cylindrical 
waves  combine  to  form  a  finite  number  of  plane  waves  propagated  in  the  dir- 
ections given  by  the  well-known  grating  diffraction  fomrula,  end  in  addition 
an  infinite  number  of  "  surface  waves  "  propagated  along  the  grating  with 
amplitudes  which  decrease  exponentially  with  the  distance  from  the  grating. 

Although  the  final  results  have  been  derived  previously  by  several 
writers'-  -' ,  our  T)rocedure,  which  is  en  extension  of  our  earlier  work  on  the 
finite  grating'--',  yields  a  direct  derivation  of  the  mode  representation  from 
the  cylindrical  wave  representation. 

II.  Perpendicular  Incidence 

We  consider  a  grating  consisting  of  an  infinite  number  of  infinitely 
long  cylinders  of  equal  circular  cross-section,  arranged  along  the  y-axis  at 
equal  intei-vals  V  ("grating  constant  "  ). 

A  plane  wave 

(1)  w^^'  =  a  exp(-iout  +  ikx), 

0 

where 

(la)  k  =  u)/c  =  2tt/\, 

impinges  perpendicularly  on  the  grating  (see  Figure  l).  The  response  of  the 
cylinder  s  (s  runs  from  -oo  to  +00)  can  be  written  in  the  following  form: 

^+00         /^  \ 

(2)  \|f   =  I_    i  A  .  H^'-^'Ctr  )  exp(in«  )  •  exp(-iujt)  . 

n=-oo 

Here  H^  '  is  the  Hankel  function  of  the  first  kind  of  order  nj  r  is  the  dis- 
n  s 

tance  of  the  field  point  P  from  the  center  of  the  cylinder  e\  Q     is  the  angle 

-9  ^ 

between  the  vector  r  and  the  positive  x-axis.  The  coefficients  A  are  determin- 

8  n 

ed  by  the  boundary  condition  at  the  surface  of  the  cylinders;  they  are  independent 


-la- 


t 


Pigore  1:  Wave  incident  on  the  infinite  grating.  Normal  incidence. 
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of  s,  since  the  grating  is  infinite.  In  papers  "by  V.  Twersky'-^-^  and  S.  Karp'--^, 
in  which  the  effects  of  multiple  scattering  are  talcen  into  account,  these  coef- 
ficients are  discussed.   In  the  present  paper  we  need  not  specify  the  A  since 
we  are  concerned  only  with  candying  out  the  sumfflation  over  s  explicitly, 
Prom  Figure  1  we  see  that 


r„  =   /x2  +  (y  - 


(3)  r„  =   Vx^  +  (y  -  sV) 


tan.,  -   ^t^ 


Where  i  and  y  are  the  coordinates  of  the  field-point  P.  We  consider  the  field 
only  for  positive  z.  The  total  scattered  wave  is 

(Zt)  y(scat)  ^  ^-l«.t  .  g^ 

where 

(4a)        S  =  1:'°°  E'°°i''\Hf^^  (k  /x2  +  (y  -  aV)^  )  .  exp(inO J  . 

8=-oo  n=-oo 

Because  the  numher  of  cylinders  in  the  grating  is  infinite,  S  must  "be  jjeriodic 
in  y,  with  the  period  V,  Hence 

(5)         S  =   H  CJt)   exp(  ^  )  . 

inr=— CO 


where 
(5a)         CU)     =  i  .  /•  ■  2  S  .  e^  (  -  ^^^  )  dy  . 


^•/. 


+  ^ 


2 


We  new  insert  the  expression  (^a)  for  S  Into  (5a);  then,  taking  into 
account  (3)  and  introducing  the  new  variehle   r^   =  y  -  sV,  we  get 


+00   +00  _      /,+  ^sY 


s=-Qon=-oo     y/    y  .^ 

*   exp  rinarctan(  ^  )  -  — y — *"  J  t 


'-^8V 


or,  equlvalently, 

+00  _+oo 


n=-oo  v 

-00 

(6) 

•  exp   (inarctaii(  J  )  -    -^I^) 

We  introduce  "^    =   ?      and  use 

arctan  S       =    it  ^°^    (i"Ht)  » 

therefore 

/  '     \^ 

exp   ( in  arctan  F    )     =    /      • ^   j  , 

and  we  obtain  finally: 

+00 

(7)  C(r)     =     ^      E      i^A    P, 


m^  y      *—  an 

n=-oo 


where 


-00  V/l+T^/ 


Here 


(7a)  q     =     kl       i  p     =       — ;^     . 


Tor      n  =  0,  one  has 


(8) 


^    expd  V q^  -  p^   )   .  for   lq|  >  jpl 


/?T  ^ 


q     -  p 

exp(  -   /p^-   q'   )    ,        for    IpI   >  |q| 


^^    ( -  >/?:: 


Yp     -  q^ 


-ij- 


(see  [5],  p.  ^7). 


In  order  to  compute  P  for  n  >  0,  we  write 


(9) 

where 

P      » 
n 

\    = 

r=0                              dp 

(10) 
that  is. 

.+00                    H^^'  (q  /l  +  !  ^     ) 

•e-^Pf     . 

(lOa) 

00                    h'I'   (q   /l  +  f  2  ) 

•  cob(p 

This  integral  can  easily  be  fotmd  from  a  general  formula  given  hy  Magntis  and 
Oberhettinger  ([5],  p.  53).  We  have 

/2^  ^  ^^^ 


0,=      1^.(^777")      '.h(^>        (/777).        ifkl>lp| 

q. 

(11)  T 

\  ■  "^  •  "=  •  (V?7) "'  •  H<^'  (1  V;v).  if  IpI  >  i< 

rrom  (9)  and  (11 )  it  follows  that 

(12)  P«  =  ^„  • 

0     0 

which  is  In  agreement  with  (8), 

Tor  n  =  1,2,3... •  computatione  based  on  (9)  and  (11 )  yield: 

(13)  ^2  -  h  ^^'^  '   ^  («?-i  >^ 

<l 


g ' 


3     J  '  cr 


for  Idl  >  IpI  , 


_<_ 


where 

For  |q|  <  IpI  we  obtain 

(13.)        P,  .  -2|e-<^'   .  <^"  (^*l)^ 

q. 


where 


I      /~"2    2 


Tte   results  (l3)  end  (13a)  may  be  generalized  Immediately  as: 
(15)        P^  =   ^  e^"^  .  cr^-^  .  (^-  i)'*  .        for  |ql  >  |p|  ; 

(15a)        P^  =  .  ^  e-  '^'  .   cr'"""^  .  (  E.^  l)n  ^      ^^^  |q,  <  |p|  ^ 

q  0^ 


As  a  matter  of  fact,  the  definition    (?)  of  P  leads  to  the  recurrence  formula 


n 


(16)        ^in  +  n  p   =  p  ^  -  ^^ 
c  q     q   a     n-1      °  p 

if  one  nses  the  relation 


(17)         H  (O  =  -  f  H  (O  +  H  -(z)  . 
n.        z  n       n— 1 


For  the  case  |q|  >  |p|,  let  us  now  assume  that  the  following  relation  holds; 


-6- 

(18)              <l'^^  •  Vl"  2.  V-"  •  (P-i*^)""^. 
with  , 

(see  eq.  (15)).  Then,  from  (l6)  and  (18)  we  have 


(19)  or 


dor 


['■^J  =  2  ^Cp-K.)"-^  •  I<r  -  E*  l(p  +  1  H.)  j 


The  relation  (19)  is  indeed  valid  if  we  put 

i<r       „ 

(20)  qV  =  2  V-  (  P-  i<^)   • 


Since  P  is  known  from  (8),  and  since  we  have  just  shown  that  (15)  is  true 

o 
for  m-l  if  it  is  true  for  n.  It  follows  "by  mathematical  induction  that  (15) 

is  true  for  all  values  of  n.  Equation  (l5a)  can  he  proven  in  the  same  way. 

In  order  to  evalus.te  S  as  given  by  (5)  we  must  consider  the  two  ceses 
Iql  ^  Ipl  separately, 

1,  Case  |ql  >  |p|.  Introducing  the  diffraction  angle  a  "by  the 
relation: 

(21)  «i^%  =  i?r    "  ""y  '  t 

(see  (7a)  )t  we  cat^  write  P  in  the  form 

(22)      T      =    S:  exp(ikx  cos  a  )  •  (cos  a  )"^^  •   (tan  a  -  i)'*  , 

n      DC  B  la  in 


-7- 


which  Is  valid  for  positive  and  negative  m,  including  m  =  0,  subject  to  the 

condition  |«i|  <  r  .  Inserting  (22)  into  (7)  and  the  resulting  expression 

for  C  (x)  into  (5) •  we  can  write  S  as  follows: 
in 


(23)  S     =     ^ 


y      ^    •    exp(ik  Jx  cos  a+  y  sin  a 


■♦■00 

n=-oo 


Reintroducing  the  time  factor  e~  "  [see  (^)J  ,  we  see  that  the  scattered  wave 
appears  as  a  strperposition  of  a  finite  numher  of  plane  waves  propagated  along 
the  directions  given  "by  the  tisual  grating  formula  (21),  The  amplitudes  of 
these  waves  are: 

?  is? 

(2^>        \  "  kVcosa    51   A  .  exp(lno^)  ^ 

m    n=-oo 

Besides  these  plane  waves  there  still  exists  an  infinite  number  of 
scattered  waves  of  a  quite  different  type.  These  correspond  to  the  case 

|p|>k|. 

2,  Case  |pl  >  |q|.  We  put 


(25)        I  =  m  7  =  ±  cosh  p^  , 


where  Im|  >  -  (  the  (+)  sign  is  taken  for  m  positive  and  the  (-)  sign  for 
m  negative).  Then  from  (5),  (7)  and  (l5e,)  we  obtain 


gj   I  np>4-oD    exp(-kx  sinh  p  )  n=»foo 

S  »  -  ^7  ]  E   .   il^^-p •  exp(iky  cosh  Pj  H   i\  exp(np^) 

I    m  >  V/\  m  n=-oo 

(26) 

UKv/X       exp(-kx  sinh  p„  ^^  .     vn  /     .   v 

^1>^/^  £     .     exp(-iky  cosh  P^)    Z.       ("^^  A^  •   e:cp(-nP^) 


■*■    >  sinh  P 

*—  m 


_  n=-oo 
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where  in  the  first  stun  over  m,  m  is  positive,  and  in  the  second,  m  is  negatire. 

— iait 
Eeintroducing  the  time  factor  e    ,  we  see  that  the  scattered  wares, 

infinite  in  n-umher,  are  !•  siirf ace  waves"  propagated  along  the  grating  in  both  y- 

directions.  Tne   spatially  damped  "  amplitudes "  of  these  surface  waves  are 


I  2i     e2p(-la:  sinh  6  )  +op 

~il50 '    H    i  \  •  exp(n3jjj),  for  positive  m 

m  n=-oo 


\    "     "  kY 


(27) 


,j              _.       eacpC-kz  sinh  p  +oo 

B^    =    -  ^  •  ^^^^ 2_  .    YZ     (-i)\-  exp(-n3^).       for  negative  m 


m       n=-oo 


....  VI  I 

in  the  directions  (+y)  and  (-y)  respectively,   Ihe  amplitudes  B   and  B  axe 


equal  if 


(28)  A_^  =  A^ 


This  is  certainly  the  case  if  the  cross-sections  of  the  cylinders  are  circular 
and  if  the  cylinders  are  perfect  conductors,  Eien  we  have 


(29)  A 


J  (ka) 


n 


n 


H^^^(ka) 
n 


(a  is  the  radius  of  the  cylinder). 

We  should  note  that  for  |m|  =  r-  "both  the  factor  cos  a  in  the  denomina- 


^  I      ti 


tor  of  B  (see  (Zk))   and  the  factor  sinh  8  in  the  denominator  of  B  and  B 

(see  (27))  hecome  zero,  hecause  in  this  case  a  =  T  »  P  =  0«  This  would  make 

t        If  m       ^        m  , 

the  amplitudes  B  ,  B   ,  B       infinite.     It  has  heen  shown  "by  IlVersky  and  EarpI 

d   m   m  ^ 

([3]  and  [}*])   that  this  "resonance  difficulty"  is  a  result  of  neglecting 
multiple  scattering. 


III.  Oblique  Incidence 

It  is  easy  to  extend  the  ahove  discussion  to  the  case  of  ohlique 
incidence  of  the  primary  wave.  Let  us  assume  (see  Figure  2)  that  the  primary 
wave  is  incident  on  the  grating  at  an  angle (-P)  with  the  x-axis.  Then 


-8a^ 


ligare   2:   Wave  incident  on  the  infinite  grating.  Oblique  incidence. 
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and  the  total  scattered  wave  'becomes 

(31)  y(scat)     ^     g-liut    ^  2 

where 

S     =    ^      JI      i'^A^e^^  •  H^-^^(k  /A(y-8Y)^)  e35)(iii0g-ilrsY  sin  p)   . 


Again,  ^     is  the  angle  which  r     makes  with  the  x-axis,  so  that 
^=       '     s  s  ' 


(32)  tan  «      = 


y  -  sV 


s  X 


In  analogy  with  (5),  we  now  have 

(33)    s  =  =  -"^  =i->  P  .  E  "  ".(-)  »^(  ^^  ) 

ia?=-oo 

with  ^ 


We  insert  S  from  (31)  into  (3^)j  then,  taking  into  account  (32),  and  introducing 
the  new  variable 


(35)  I  =  21^ 
we  oTstain  finally 

+00  .     - 

(36)  C^(x)    =    f   E      i\  e^^P     .    P^ 

n=-oo  * 

where  /  v 

(37)     p,  =  /"df  .  =^1)  6  /r7?)(LiiL| " .  ,-ip., 


-10- 


Here, 


a    = 


i     p       =     I 


[ 


2ti!ii 


-  k  sin  P 


1 


(38)    <  that  is. 


p      a    p  -  q  sin  P    =     q.  (  ^  -  sin  P  )       , 


where  p  is  defined  by  (7a). 

The  only  difference  "between  (7)  and  (37)  is  that  p  has  heen  replaced 
t 
"by  p  •     Hence  from  (15)  we  may  conclude  immediately  that 


(39) 


?n    =    7    ^"^    •    ^''"'-     ^#-^^' 


J         _  1^-1  • 


for  \q\  >  |p  1      , 


^   ^n    =    -T^""^    •     ^  •    (|t*1)''  forU|<|p'|      . 


where 
(39a) 


5-  =  /777 ;      ^'  =   / 


'2        2 

P      -  <1 


Jgain,  in  order  to  evaluate  S  we  consider  the  two  cases   jqj  ^  jp   |   separately, 
1,  Case   |q|  >  Ip   1.     We  put,  in  analogy  with  (21): 


(ko) 


^  =     ^    -  sin  p     =     sin  o^     . 


so  that 
(40a) 


sin  a      +     sin  B 


=     m 


V    • 


Then  jq]  >  |p   |   implies  that 


(41) 


1    >     1^    -     Sin  Pi, 


that  is. 


1   sin  a\ 


sin  a^l     <    1 
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From  (36]_j  (39)  and  (33)  we  find 

(42)  S     =    ^JI    ~  •  expjik  Jx  cos  a+y  sin  a  llJI    A^«  exp[in(a+p)]   . 

m   cos  a      V   I  jyn=-OD 

The  values  of  o  range  from  m  ^  -  j-  (1  -  sin  g)  to  m  _  r  (1  "••  sin  P). 

Reintroducing  the  time-factor  e""   we  see  that  the  scattered  wave  appears 

as  a  superi^osition  of  a  finite  n\mber  of  plane  waves  propagated  eQong  the  directions 

_  *  II 

given  by  the  usual  grating  fromula  (40a).  Die  angle  a  is  the  "diffraction  angle  , 

that  is^  the  angle  between  the  z-axis  (the  normal  to  the  grating)  and  the  propaga- 
tion vector  of  the  diffracted  wsvea.  The  aaiplitudes  of  these  waves  are: 

_  ^  +00 

(^3) 


(W+) 


B 
m 

2             .     5^ 

^  '^'^  \          too' 

2. 

t, 
Case     q    <    p    j.     Put 

^=     m^    -     sin  p 

—      m 


where  the   (+)  sign  is  taken  for  m  >  r-  (1+  sin  P),  and  the   (-)  sign  for 
m  <  -  ^  (1  -  sin  p).     Then  from   (33),    (36)  and   (39)  we   obtain 


{ 


-       ,  exp(-kx  sinh  "p  ) 

S     =    -  i    ]  r  •  ——Z ^  •  expdl^  cosh  PJ 

i!^(l+9in  P)  ^^  ^ 


+  00         _  -  _     _  exp(-!ac  sinh  P   ) 


(45)  •  E        i\  exp[(ip  +  p^)n]  +  J2 


m 


nK  V/\(l-8in  p)  ^m 


+  00  _      n  1 

.   exp(iky  cosh  ?  )  J^         (-i)^  .   A    e3cpr(ip  -  p  )n|    I 

"  n=-oo         ^  m  J  j' 


This  is  the  generalization  of  (26)  for  oblique  incidence  of  the  primary  wave. 
Reintroducing  the  time  factor  e"  *"  ,  we  see  that  the  scattered. waves. 
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infinite  in  ntupiber,  are  "  siorf ace  wares  "  propagated  along  the  grating  in  both 
y-directions.  The  mode  characterized  "by  m  propagates  with  the  velocity 


ih6)  V„  = 


m      I   X 


Im  ^  -  sin  Pl 


where  the  denominator  is  greater  than  1.  Its  wave  length  is 


(^7)  A 


"^      1  m^-  sin  Pl 


and  its  arnplitnde 

—r  _.        exp(-Ia:  sinh  ?„)         +op  .        _      - 

B_     =     -^'  Z ^*    II    iX  ®^L(^^PmH»         for  positive  m, 


(^48) 


sinh  P  n=-c» 

IQ 


„         21    «^^-^  ^^^  K^     i3> 


B^    =     -  ^  •  :: 2^  •    H     (-i)''.A^e2p[(ip-p^)n],  for  negative  m. 


sinh  P  n=-oo 
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